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Abstract 

We consider the double sinh-Gordon potential which is a quasi-exactly solv- 
able problem and show that in this case one has two sets of Bender-Dunne 
orthogonal polynomials . We study in some detail the various properties of 
these polynomials and the corresponding quotient polynomials. In particular, 
we show that the weight functions for these polynomials are not always posi- 
tive. We also study the orthogonal polynomials of the double sine-Gordon po- 
tential which is related to the double sinh-Gordon case by an anti-isospectral 
transformation. Finally we discover a new quasi-exactly solvable problem by 
making use of the anti-isospectral transformation. 
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I. INTRODUCTION 



Recently, in an interesting paper, Krajewska et al. [|T| have introduced an anti-isospectral 
transformation ( which they also called as duality transformation ). In particular using this 
transformation they relate the spectra of quasi-exactly solvable ( QES ) potentials Vi and 
V2 given by 

Vi{x) = x^{ax^ + hf - na{2M + ?,)x^ (1.1) 
V2{x) = x\ax^ - bf - ha{2M + 3)x^ (1.2) 

where a, 6 > and M is nonnegative integer. It may be noted here that the duality property 
does not hold for other (exactly non-calculable) levels. It is clearly of great interest to 
understand this new transformation in some detail and explore it's various consequences. 
For example, one would like to know whether one can discover new QES problems by using 
this transformation. Secondly are the QES levels of the dual potentials related even if they 
are valid over different domain? Further, whether the number of QES levels in the two cases 
are identical or not. 

Another recent development is the work of Bender-Dunne |Q where they have shown that 
the eigenfunctions of the Schrodinger equation for a quasi-exactly solvable ( QES ) problem 
is the generating function for a set of orthogonal polynomials {P„(i?)} in the energy variable 
E. It was further shown in one specific example that these polynomials satisfy the three-term 
recursion relation 

P^{E) = EPn-,{E) + CnPn-2{E) (1.3) 

where C„ is E independent quantity. Using the well known theorem " the necessary 

and sufficient condition for a family of polynomials {Pn} ( with degree Pn = n ) to form 
an orthogonal polynomial system is that {Pn} satisfy a three-term recursion relation of the 
form 

P^{E) = [Ar^E + B,,) Pn-l{E) + CnPn-2iE), U > 1 (1.4) 
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where the coefficients An, Bn and C„ are independent of E, An 7^ 0, Ci = 0, C„ 7^ 
for n > 1", it then followed that {Pn{E)} for this problem forms an orthogonal set of 
polynomials with respect to some weight function, w{E). Recently several authors have 
studied the Bender-Dunne polynomials in detail In fact it has been claimed that 

the Bender-Dunne construction is quite universal and valid for any quasi-exactly solvable 
model in both one as well as multi-dimensions. 

However, in a recent note we have discussed three QES problems for all of which 
the Schroodinger equation can be transformed to Heun's equation and further in all these 
cases the three-term recursion relation satisfied by the Bender-Dunne polynomials is not 



of the form as given by Eq. (|1.4| ) and hence does not form an orthogonal set. It is worth 
pointing out that in all these cases the Hamiltonian can not be written in terms of quadratic 
generator of the Sl{2) algebra. We suspect that this may be the reason why the polynomials 
in these examples do not form an orthogonal set. 

The question which we would like to raise and study in this paper is regarding the prop- 
erties of the Bender-Dunne polynomials in case they form an orthogonal set. For example in 
the example discussed by Bender-Dunne (as well as in most other QES examples discussed 
so far), for a given potential, the QES states have either even or odd number of nodes. The 
question is what happens in a QES problem if one can obtain states with even as well as 
odd number of nodes. In particular, do the polynomials corresponding to the even and odd 
number of nodes together form an orthogonal set or if the polynomials corresponding to even 
number of nodes form one orthogonal set and polynomials with odd number of nodes form 
a separate orthogonal set ? Secondly, how universal are the properties of Bender-Dunne or- 
thogonal polynomials in such cases? For example, are the weight functions always positive? 
Do the moments of weight function have pure power growth? |Q 

The purpose of this note is to explore in some detail the various issues raised above. 
In particular, in Sec. II, we discuss the double sinh-Gordon (DSHG) model which is an 
QES problem but where states with odd as well as even number of nodes are known for 
a given potential. We show that in this case there are two sets of orthogonal polynomials 



with polynomials corresponding to even number of nodes forming one orthogonal set while 
the polynomials corresponding to odd number of nodes forming another orthogonal set. We 
study the properties of these polynomials as well as the corresponding quotient polynomials 
in detail and show that in many respect they are similar to the Bender-Dunne polynomials. 
However, unlike the Bender-Dunne case, the weight functions are not always positive for 
either of the orthogonal polynomial sets. In Sec. Ill, we discuss several consequences of the 
duality symmetry in the context of the QES problems. In particular, we study two QES 
problems, in both of which the domain of validity of the potential and it's dual are different. 
As a result we show that in the case of the double sine-Gordon ( DSG ) potential, the number 
of QES levels are almost half of those in the double sinh-Gordon case ( which is dual to 
it ). However, even then many of the predictions about their corresponding Bender-Dunne 
polynomials still go through. For example, the QES levels of DSG are still related to the 
corresponding levels of DSGH case. Further, the weight functions here are also not always 
positive. Finally using duality, we discover a new QES problem which has so far not been 
discussed in the literature. Sec. IV is reserved for discussions. 



II. DSHG SYSTEM AND ORTHOGONAL POLYNOMIALS 

The DSHG system is one of the few double well problems in quantum mechanics which 
is quasi-exactly solvable This double well system has found application in several 

different branches of physics starting from the theory of diffusion in bistable field , 
quantum theory of instantons 0,|12|, quantum theory of molecules and also as a model for 
nonlinear coherent structure |]13|,|1^. The Hamiltonian for this case is given by ( we shall 
assume h = 2m = 1 throughout this paper) 

(P 

H = -— + (C cosh 2x - Mf (2.1) 

where C is a positive parameter. The value of M is not restricted in principle but it has 
been shown |^ that the solutions for first M levels are exactly known in case M is a positive 
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integer. Note further that for M > C, this potential has a double well structure with the 
two minima lying at cosh2xo = Let us now derive the three-term recursion relation in 



this case. On substituting 



(2.2) 



in the Schroodinger equation Hi/j = Eip with H as given by Eq. we obtain 



"(x) - 2C sinh 2a;(^'(x) + {E - ~ C ) + 2(M - 1)C cosh 2x (/)(x) = 



(2.3) 



On further substituting 



z = cosh2x — 1; (p = z'^ 



Rr,{E) fz + 2 



n=0 



(2.4) 



we obtain the three-term recursion relation {n > 0). 



Rn+2{E) - [n^ + 4(s + On + 4s' + E - - - 2(M - l)c] Rn{E) 

- AC[M +l-2s-n]n{n-l)Rn-2{E) = (2.5) 

provided 2s' = s i.e. either s = or s = |. Thus we have two sets of independent solutions 
: the even states (i.e. states with even number of nodes) for s = and the odd states for 
s = \. Note that unlike the Bender-Dunne (or most other QES) cases, s is not contained 
in the potential and this is perhaps related to the fact that for any integer M (> 2) the 
QES solutions corresponding to both even and odd states are obtained. From Eq. ( |2.5| ) we 
observe that the even and odd polynomials Rn{E) do not mix with each other and hence 
we have two separate three-term recursion relations depending on whether n is odd or even. 
In particular, it is easily shown that the three-term recursion relations corresponding to the 
even and odd n cases respectively are given by (tt, > 1) 



Pn{E) 



Qn{E) - 



An^ + 8n(s + C - 1) + 4s' - 8s + 4 - 6C + £^ - (M + C)^] Pn~i{E) 
8C(n - 1) (2n - 3) (M + 3 - 2s - 2n) Pn-2 = (2.6) 

4n' + 4n(2s + 2C - 1) + 4s' - 4s + 1 - 2C + ^ - (M + C)'] Qn-i{E) 
8C{n - l){2n - 1) (M + 2 - 2s - 2n) Qn-2 = (2.7) 



and 



with Po{E) — 1, Qo{E) — 1. These recursion relations generate a set of monic polynomials, 
of which the first few are 

Po{E) = 1 

Pi(E) =£ + 4s2 + 2C 
P2{E) ^£^ + £ [Ss^ + 8s + 12C + 4 

+ (45^ + 4s + IOC + 4) (45^ + 2C) + 8C (M - 2s - 1) (2.8) 

QoiE) = 1 

Qi(E) =£ + 4s2 + 4s + 1 + 6C 
Q2{E) =£^ + £ [8s^ + 8s + 20C + lO' 

+ (4s2 + 4s + 6C + 1) (4s2 + 4s + 14C + 9) + 24C (M - 2s - 2) (2.9) 

where 

£ = E-{M + Cy (2.10) 

Note that unlike the Bender-Dunne case, the polynomails Pn{E) and Qn{E) are not eigen- 
functions of parity. It is easily seen that when M is a positive integer, exact solutions for first 
M levels are obtained. In particular, if M is odd (even) integer, then solutions with even 
number of nodes (s = 0) are obtained when the coefficient of P„_2 iQn-2) vanishes. Similarly 
if M is odd (even) integer, the solutions with odd number of nodes (s = |) are obtained 
when the coefficient of Qn-2 (-fn-2) vanishes. Further for M even (say 2A;-|-2, A; = 0, 1, 2 • • • 
), half the levels, i.e. {k+ 1) levels are obtained each from the zeros of the orthogonal poly- 
nomials Pk+i{E) and Qk+i{E). On the other hand when M is odd (say 2A;-|-1, k = 0,1 ■ ■ ■) 
then k + 1 and k levels are obtained from the zeros of the orthogonal polynomials Pk+i{E) 
and Qk{E) respectively. 

We have studied the properties of the polynomial sets {P„(£')} and { Qn{E)} and we 
find that many of their properties are almost the same and in turn they are similar to the 

6 



Bender-Dunne polynomials . In particular, since in both the cases the recursion relations 
are similar to those given by Eq. ( p. .41 ) , hence for all values of the parameters M and s, 
they form an orthogonal set . Secondly, the wave function iplx, E) is the generating function 
for the polynomials Pn{E) as well as Qn{E). Thirdly when M is a positive integer, both of 
these polynomials exhibit factorization property whose precise form depends on whether M 
is even or odd. In particular, ifM = 2A; + l, fc = 0,l,2-- - then one has {n > 0) 

Pk+l+niE) = Pk+l{E)Pn{E) 

Qk+n{E) = Qk{E)Q^{E) (2.11) 
On the other hand if M = 2fc + 2, A; = 0, 1, 2 ■ ■ -, then 

Pk+l+n{E) = Pk+l{E)Rn{E) 

Qk+i+n{E) = Qk+i{E)Sn{E) (2.12) 

where Po(^) = <5o(^) = ^o(^) = Sq{E) = I. Following Ref. § it is easily shown that the 
polynomial sets {Pn{E)}, {Qn{E)} , {_R„(£')} and {Sn{E)} correspond to the non-exact 
spectrum for this problem. 

To illustrate this factorization we list in factored form, the first few polynomials of both 
the types in case M = 3 and M = 4. 

M = 3 



P,{E) = 1 
Pi{E)=S + 2C 

P2{E) =£^ + £ [12C + 4] + 20C' + 24C 
P3(^) = [^ + 18C + 16]P2(^) 

Pi{E) = + £{A6C + 52) + (36 + 28C)(16 + 18C) - 240(1 P2iE). 



(2.13) 



Qo{E) = 1 
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Qi{E)=£ + QC + A 



Qi{E) 



Q^{E) = [S^ + £{36C + 52) + (36 + 22C)(16 + 14C) - 160q Qi(^). (2.14) 

where S is given by Eq. ( |2.(j| ). Let us notice that in this case P2{E), is a common factor of 
Pn{E) for n > 2. The zeros of P2{E) are at 



Eo = 7 + C-2Jl + AC 



E2 = 7 + C + 2V1 + 4C 



(2.15) 



which give the ground and second excited state eigenvalues for the potential when M = 3. 
On the other hand Qi{E) is a common factor of Qn{E) for n > 1 . The zeros of Qi{E) are 
at 



E, = e + 5 



(2.16) 



which give the first excited state eigenvalue for the same potential (i.e. when M = 3). 
In this way for M = 3 one obtains exact energy eigenvalues for the first three levels. 



M = 4 



PoiE) = 1 
Pi{E)=S + 2C + l 

P2{E) =S'^ + S [12C + 10] + 20C^ + 44C + 9 
P3(^) = [^ + 18C + 25]P2(^) 

Pi{E) = + ^(46C + 64) + (28C + 49)(18C + 25) - 240(1 P2{E). 



(2.17) 



QoiE) = 1 

Qi(E) =^ + 6C + l 

Q2{E) =S^ + S [20C + 10] + 84C' + 116C + 9 
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Q^{E) = [£ + 22C + 25] Q2{E) 

QiiE) = [S^ + S{52C + 74) + (30C + 49)(22C + 25) - 336C] Q2{E). (2.18) 
In this case P2 {Q2) is a common factor of Pn (Qn) for n >2 . The zeros of P2{E) are at 



Ei = C' + 2C + ll-4Ve + C+l 



i?3 = C + 2C + 11 + + C + 1 (2. 19) 

which give the first and the third excited state energies for the potential when M = 4. On 
the other hand the zeros of Q2{E) are at 



i?o = c'-2c + ii- Vc'-c + i 

£^2 = C' - 2C + 11 + - C + 1 (2.20) 

which give the ground and the second excited state energies for the same potential (M = 4). 
The corresponding eigenf unctions can be easily obtained by evaluating iIj^x) in Eq. ( p.4|) 
at these values of E and not surprisingly they are the same as given by Razavy p| ( with 
appropriate change of parameters ). 

The norms (squared) of the orthogonal set of polynomials Pn{E) and Qn{E) can be easily 
determined by using the recursion relations and we find ( assuming 70^'' = 70'^'' = 1) 

n 

7W = (-8C)"ri! n [(2A; - 1) (M - 2s - 2fc + 1)] 

k=l 
n 

7(9) = (-8C)"n! n [(2A; + 1) (M - 2s - 2k)] (2.21) 

k=l 

We observe from here that 7^^-' vanishes for 2n > M — 2s + 1 while vanishes for 2n > 
M — 2s provided M is a positive integer and as remarked by Bender-Dunne, this vanishing of 
norm is an alternative characterization of the QES problem. It may be noted that unlike the 
Bender-Dunne case, the norms are alternative in sign for 2n < M — 2s + 1 or 2n < M — 2s 
depending on the set {/-'„(£')} or {Qn{E)} respectively. 

Using the factorization property as given by Eqs. ( |2.11| ) and ( |2.12| ) and using the recur- 



sion relations ( P^ ) and ( p^.T] ) satisfied by P„ and Qn, it is easily shown that F„ Q„ i?„ and 
Sn satisfy the following recursion relations 



Pn{E) = [(M + 1 + 2nY + 4(M + 1 + 2n)(s + C - 1) + - 8s + 4 - 6C + ^] Pn-i{E) 

+ 4C(M - 1 + 2n)(M + 2n - 2) [2 - 2s - 2n] Pn-2{E) 
Qn{E) =[{M-1 + 2nf + 2(M - 1 + 2n)(2s + 2C - 1) + 45^ - 4s + 1 - 2C + ^] Qn-i(^) 

+ 4C(M - 3 + 2n)(M + 2n - 2) [3 - 2s - 2n] Qn-2{E) 
Rn{E) = [(M + 2n)2 + 4(M + 2n)(s + C - 1) + 4s2 - 8s + 4 - 6C + S] Rn-i{E) 

+ 4C(M - 1 + 2n)(M + 2n - 3) [3 - 2s - 2n] Rn-2{E) 
SniE) = [{M + 2nf + 2(M + 2n)(2s + 2C - 1) + 4s^ - 4s + 1 - 2C + ^] Sn-i{E) 

+ 4C(M - 2 + 2n)(M + 2n - 1) [2 - 2s - 2n] ^„_2(^) (2.22) 



Using these recursion relations, it is straightforward to obtain the norms of theses poly- 
nomials and show that they are all positive. In particular, we obtain 

n 

In = (40" n [(M + 2k + l){M + 2k){2k + 2s)] 
fc=i 

7? = (40" n [{M + 2k){M + 2k- l){2k + 2s - 1)] 

k=l 
n 

In = (40" n [(^ + 2k + l){M + 2k- l){2k + 2.3- 1)] 

k=l 



tI = (40" n [(^ + 2k + l){M + 2k){2k + 2s)] 

k=l 



(2.23) 



Note that in Eqs. (|2.22| ) and ( 2.23|) M and s take only specific values. For example, in the 
case of Pn{E) M is odd and s = 0, while in the case of Qn{E) M is odd and s = |. On 
the other hand, in the case of Rn{E) M is even and s = | while in the case of Sn{E) M 
is even but s = 0. 

Finally, we can also obtain the weight function w{E) for our polynomials {Pn{E)} and 
{Qn{E)} by using the expression as derived by Krajewska a/ j^. In particular, the weight 
function w^^\E) can be written as 

M 



w^P\E) = J2w^^>5{E-Ek) 

k=l 



(2.24) 
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where the numbers w^^^ satisfy the algebraic equation 

M 

Y,PniEk)w'i'^ = 5nO (2.25) 

k=l 

A similar relation also exists for the weight function of the Q polynomials. 

As an illustration we have computed the weight functions ly^^-* and w^^ in case M = 3 
and 4. For example, when M = 3, using Eqs. (^I^ ) to (g]^) and ( |ZT^ ) we find that 

(P) _ 1 2C + 1 







2 2 VI + 4C^ 
(p) _ 1 2C + 1 

"2 + 2v/r+4e 

= 1 (2.26) 



On the other hand, when M = 4, using Eqs. (|]2D, (|]2|),(|2]T3) and ( ^IsD we find 



(.) 1 c 

U^n - 



W 



2 2VC' - C + 1 

(p) _ 1 _ C + 1 

2 2VC + C+1 



^2 = o H 



2 2VC' -C+1 

As a cross check on our calculation of weight functions, we have calculated the norm of the 
polynomials by using the basic relations 



In 



J dEw^^\E)[P.^{E)f- 71'?) = j dEw^'^\E)[Q^{E)f (2.28) 



and for M = 3 as well as 4 we have verified that we get the same answer as given by Eq. 
( p. 21 ). It is worth pointing out that unlike the Bender-Dunne example, in none of our cases 



the weight functions are always positive. Actually this is not all that surprising. As has 
been proved by Finkel et al. , if the three-term recursion relation is of the form 

Pk+i = {E-h)Pk-akPk-i, k>0 (2.29) 
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with ao = and a„+i = 0, then the weight functions are all positive if bk is real for < A; < n 
and ajt > for 1 < k < n. On comparing Eq. ( |2.29D with our recursion relation ( p.6| ) and 



( p.TD we find that in both of our cases, ak is in fact < for 1 < k < n. 

Once the weight function is known, then one can also calculate the moments of w{E), 
defined by 

/i„ = y dEw{E)E'^ (2.30) 

Since in our case, the polynomials {P„(i?)} and {Qn{E)} are not eigenfunctions of parity, 
hence unlike the Bender-Dunne case, odd moments do not vanish in our case. For example, 
the first few moments in our case for M = 3 are ( we choose = yUp'^'' = 1 without any 
loss of generality ) 

^^^r^ = (3 + C)^ - 2C; /i?^ = -16C + [(3 + C)^ - 2Cf 

= (3 + - 6C - 4; /i^'^) = [(3 + 0' - 6C - 4]^ (2.31) 

Following Finkel et al. , we can also calculate the growth rate of the moments. In partic- 
ular, since in both of our cases, 7^ for 1 < A; < n where ak as given by Eqs. ( |2.29| ), ( ^.6| ) 
and P77| , hence following their discussion it is easily shown that for large n, to leading order 

/z(f), /i(f)~(M + C)'" (2.32) 

Thus in our case the moments have a pure power growth. 

III. NEW QES POTENTIALS FROM ANTI-ISOSPECTRAL 
TRANSFORMATIONS 

In an interesting paper, Krajewska et al. |^ have recently discussed the consequences of 
the anti-isospectral transformation ( also termed as duality transformation ) in the context 
of the QES problems. In particular, they have shown that under the transformations x — *■ 
ix = ?/, if a potential V{x) goes to Vijj), i.e. 
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V{x) V{y) (3.1) 

then the QES levels of the two are also related. In particular, they have shown that if 
M levels of the potential V[x) are QES levels with energy eigenvalues and eigenfunctions 
Ekik = 0, 1, ■ ■ ■ , M — 1) and -0^(3^) respectively then the energy eigenfunctions of V{y) are 
given by 

Ek = -EM-i-k, '>Pk{y) = i^M-i-kiix) (3.2) 

As an illustration, these authors have discussed the x^- potentials as given by Eqs. ( [L.l| ) 
and ( pr^ ) and explained in details how the eigen spectra of these two dual potentials are 
related. However in this particular case, the domain of x and y are the same. 

The purpose of this section is to explore the consequences of this symmetry when the 
domain of validity of V{x) and V{y) are different and to see how many of the results derived 
in Ref. go through in this case. In particular, do V{x) and V{y), hold the same number 
of bound states ? Are the bound states of the two potentials still related by Eq. (|3.2| )? To 
begin with, we would also like to comment that, the duality symmetry is only useful in case 
the potential V{x) is symmetric in x as otherwise the dual potential V{y) will be a complex 
potential . 

A. DSG Example 

Let us consider the potential corresponding to the DSHG case, and try to explore the 
consequences of the duality symmetry in this case. On applying the duality transformation 
X ^ i9 to the Schroodinger equation Hi/j = Eip with H as given by Eq. ( p.l|) we obtain the 
following Schrodinger equation for the DSG equation. 

■^{9) = Eilj{e), E = -E. (3.3) 

Note that whereas the domain of x in DSHG equation is — oo < x < oo in the DSG case, 
the domain is < ^ < 27r provided we are solving the problem on a circle as we do here. As 
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- {Ccos29-Mf 



d9^ 



a result, unlike the DSHG case, in the DSG case due to the periodicity constraint, one has 
to demand that the wave function tp{9) must be invariant under 6 ^ 9 + n i.e. 

^l;{e + tt) = ^{9) (3.4) 

Looking at the QES eigenfunctions of the DSHG case, it is easy to see that when M is even 
{i.e. M = 2,4- ■ ■) then the eigenfunctions of the DSG case do not remain invariant under 
9 —>■ 9 + 71 but they change sign. As a result, the QES levels of DSG when M is even must 
be rejected. Thus on physical ground, we find that in case the domain of V{x) and V{y) 
are different then the number of QES levels of V{x) and V{y) need not be identical. In 
particular, boundary conditions may forbid some wave functions to be eigenfunctions. For 
example, for M = 2, naively one would have two exact energy eigenstates of DSG equation 
as given by 

ijj^^^{9) = sinee-2^°^2e 



However both of these are unacceptable eigenfunctions of DSG, being odd under 9 ^ 9 + 7r. 
Thus we have shown that even though DSG and DSHG are dual system, DSG equation 
has approximately only half the number of QES levels compared to the DSHG case. In 
particular, when M is an odd integer (M = 1, 3, 5 ■ ■ ■) only then the DSG is a QES system 
and in that case the first M levels are exactly known, and in fact they can be immediately 
obtained from the corresponding DSHG case by making use of the duality relation ( p. 21) . 
For example, for M = 1 the ground state of the above DSG is given by 

^o = -(l + C') ^o = e-^™^2'^ (3.6) 
On the other hand for M = 3, the first 3 levels of the DSG equation are are given by 



g- 1 cos 29 (3 7^) 



Eo = -7-C-2^1 + 4e; 
^i = -C'-5; ^/^i = sin2^e-2=°^2e 



2C- 1 VI +4(2- U cos 2^ 



^2 = -7-C' + 2Vl + 4C2; 



2C+ |v^lT4C2 + l|cos20 



14 



The energy eigenstates of the DSG equation have recently been obtained by Habib et al. 
and not surprisingly our results are the same as obtained by them. 

One can now study the Bender-Dunne polynomials of the DSG equation and it is easy 
to see that most of the discussion of the last section ( for the DSHG case ) also goes through 
in this case except that now only odd values of M give us a QES system. For example, 
the recursion relations for P„ and Qn as given in the last section are also valid in this case 
provided we replace E by ~E. Thus the three-term recursion relation for P„ as given by Eq. 
( p.6|) with E changed to —E will give an exact solution for the DSG case only if s = while 
the three-term recursion relation for Qn as given by Eq. ( p.7| ) ( with E changed to —E ) will 
give an exact solution only if s = |. Most of the other properties in the last section continue 
to be valid in the DSG case also ( with the obvious replacement of E by —E ) except now 
only Pn{E) and Qn{E) exist while Rn{E) and Sn{E) as given by Eq. (|2.22|) do not exist 
( Note that Pn and Qn do not give QES solution if M is even integer). In particular, we 
want to emphasize that the the norms of the polynomial sets {P„(E')}, {Qn{E)}, {Pn{E)} 
and {Qn{E)} are unchanged from those of the DSHG case. However the weight functions 
get interchanged . For example for M = 3 the weight functions for the DSG case are 

(P) 1 2C + 1 
Wn = t: + 



w 



2 2V1 + 4C^ 
(p) = i _ 2C + 1 
^ 2 2V1 + 4C^ 

w[''> = 1 (3.8) 

Thues even in this case, the weight function are not always positive. Finally it is easily 
shown that the moments for the DSG and DSHG cases are related by 

^^'')|D5G=(-l)>^'')b5/.G (3.9) 

SO that for large n, to leading order, the DSG moments also have a pure power growth rate 
as given by 

//„~(-l)"(M + C)'". (3.10) 
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B. New QES Potential 



We shall now point out a non-trivial application of the duality symmetry. In particular, 
using this, we obtain an entirely new QES potential which has so far not been discussed in 
the literature |1T3[ . 

Consider the following potential 



V{x) 





8 sinh^ 1) sinh' f + 2 






8 


1 + ^ + sinh^ f 


2 



(3.11) 



which is obtained in the context of the stability analysis of the (p^ -kink solution in 
1+1 dimensions [|T^. As has been shown before, in this case the Schroodinger equation 
can be converted to Heun's equation with four regular singular points. Further, for any e, 
the ground stat of the Schroodinger equation corresponding to this potential is given by 
(n = 2m = 1) 



^0 = A^o 



e^ + l 



e^ + l + e^sinh'f J ^ 



e2 sinh^ if 



+ 



(3.12) 



e^ + l + e^sinh^f ' 2_ 

Further for the special case of, = |, the second excited state is also analytically known 
and given by 



1p2 = N2 



3 + sinh^ f 



sinh^ f 



E2 = 



(3.13) 



3 + sinh' f 4_ 

Thus for = I this is a QES system. It may be noted that even though it is a QES system 
in one dimension, the Hamiltonian can not be written in this case in terms of the quadratic 
generators of S'/(2) |16|. Further, as shown by us recently 0, in this case the Bender-Dunne 
polynomials do not form an orthogonal set . We shall now apply the duality transformation 
and obtain a new QES system which has not been discussed before. 

On considering the duality transformation x ^ iO in the Schroodinger equation Hifj = 
Eip corresponding to the potential ( p.ll|) , we find that we have a new periodic potential 



V{9) 





8 sin^ f + 4(^ - 1) sin^ f + 2 






8 


1 + i - sin^ f ; 


to 



(3.14) 
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Notice that the domain of V{x) and V{9) are very different i.e. whereas — oo < fix < 



oo, < fi9 < 2n. However, the predictions of duahty symmetry as given by Eq. ( |3.2| ) are 
still valid. In particular, using Eq. ( |3.2| ) we predict that the potential V{9) as given above 
is a QES system in case = | and it's ground and second excited states energy eigenvalues 
are given by 



0, ^2 



3 - sin2 f 



+ 



sm 



2 ^ 



4 3 - sin^ f 



e2 + 1 - e2 sin^ f 



e^sin^f 



2 ■ 2 ud 

e'^ sm ^ 



(3.15) 



In particular, we would like to emphasize that for the QES potential ( p.l4|) , the second 
excited state -02 is known at all real values of e, while ipQ is only known at = |. One 
can solve the Schroodinger equation explicitly for the potential (|3.14|) and check that indeed 
these are the energy eigenstates of the system. In particular, note that both ipQ and -02 
are invariant under fiO ^ fi9 + 2t\ i.e. + 27r) = il){ii9). It is really remarkable that 

using duality we are able to obtain a new QES system. In this case one can also obtain the 
three-term recursion relation for the Bender-Dunne polynomials and show that they do 
not form an orthogonal set. This is related to the fact that the corresponding Schroodinger 
equation can be converted to Heun's equation with four regular singular points. 

It may be worthwhile to look at all known QES systems and see if the corresponding dual 
systems have already been discussed in the literature or not. In this way one may discover 
some more QES systems. 



IV. DISCUSSIONS 

In this paper we have analyzed in some detail a QES system for which the energy 
eigenstates for levels with odd as well as even number of nodes are known for a given 
potential. We have seen that in this case one obtains two independent sets of orthogonal 
polynomials. Further, we have also seen that the weight functions in this case are not 
necessarily positive. It will be interesting to study few other examples of the same type 
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and enquire if in those cases too one obtains two sets of orthogonal polynomials or not. 
Further, whether the weight functions in such cases are always positive or not. We have 
also compared the spectrum of two dual potentials in case their domain of validity are quite 
different. We have seen in one case that because of different boundary conditions, the QES 
spectrum of the dual potentials is not the same. It will be interesting to study several such 
QES dual potentials and see if some general conclusions can be obtained in these cases. 

Finally, it would be really interesting if one can discover some new QES potentials by 
making use of the duality symmetry. 



18 



[1] 

[2] 



[4] 
[5] 
[6] 
[7] 



[9 
[10 

[11 
[12, 
[13' 
[14 
[15 



REFERENCES 

A. Krajewska, A. Ushveridze and Z. Walczak, Mod. Phys. Lett. A 12 (1997) 1225. 
C. M. Bender and G. V. Dunne, J. Math. Phys. 37 (1996) 6. 

T.S. Chihara, An Introduction to Orthogonal Polynomials , Gordon and Breach, New 
York, (1975); A. Erdelyi, W. Magnus, F. Oberhettinger and F. G. Tricomi, Higher 
Transcendental Functions, Vol.11, McGraw-Hill , New York, (1953). 

F. Finkel, A. Gonzalez-Lopez And M. A. Rodriguez, J. Math. Phys. 37 (1996) 3954. 

A. Krajewska, A. Ushveridze and Z. Walczak, Mod. Phys. Lett. A 12 (1997) 1131. 

C. M. Bender, G. V. Dunne and M. Moshe, Phys. Rev. A 55 (1997) 2625. 

A. Khare and B. P. Mandal, [physics/ 9709"043i IP-BBSR/97-36. 



M. Razavy, Am. J. Phys. 48 (1980) 285. 

S. Habib, A. Khare, and A. Saxena, (In preparation). 

N. G. van Kampen, J. Stat. Phys. 17 (1977) 71. 

A. M. Polyakov, Nucl. Phys. B 120 (1977) 429. 

N. H. Christ and T.D. Lee Phys. Rev. D 12 (1975) 1606. 

A. Khare, S. Habib and A. Saxena, Phys. Rev. Lett. 79 (1997), Nov. 10 issue. 

S. N. Behera and A. Khare, J. de Phys. C6 (1981) 314. 

A. Ushveridze, Quasi- Exactly Solvable Models in Quantum Mechanics, Inst, of Physics 
Pubhshing, Bristol, (1994). 



[16] D. P. Jatkar, C. Nagaraja Kumar and A. Khare, Phys. Lett. A 142 (1989) 200. 



19 



